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Abstract 

^ ■ We introduce new symmetry-based methods to test for isotropy 

. in cosmic microwave background radiation. Each angular multipole 

QQ I is factored into unique products of power eigenvectors, related mul- 

CN ' tipoles and singular values that provide 2 new rotationally invariant 

00 ' measures mode by mode. The power entropy and directional entropy 

! are new tests of randomness that are independent of the usual CMB 

I power. Simulated galactic plane contamination is readily identified, 

^ ■ and the new procedures mesh perfectly with linear transformations 

• i-H . employed for windowed-sky analysis. The ILC -WMAP data maps 

^ ' show 7 axes well aligned with one another and the direction Virgo. 

■ Parameter free statistics find 12 independent cases of extraordinary 
axial alignment, low power entropy, or both having 5% probability or 
lower in an isotropic distribution. Isotropy of the ILC maps is ruled 
out to confidence levels of better than 99.9%, whether or not coinci- 
dences with other puzzles coming from the Virgo axis are included. 
Our work shows that anisotropy is not confined to the low 1 region, 
but extends over a much larger / range. 

Keywords: cosmic microwave background, methods:data analysis, meth- 
ods:statistical 
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1 Introduction 



Studies of the cosmic microwave background (CMB) were long posed in terms 
of the temperature power spectrum. The power spectrum is invariant under 
rotations, and by itself cannot in principle test the basic postulate that the 
radiation field should be statistically isotropic. The availability of high qual- 
ity data from WMAP (Hinshaw et al 2003, 2007) has allowed isotropy to be 
examined critically, and with surprising outcomes. 

A long-standing question exists in the unexpected size of low-Z multipoles. 
Interpretation of this is stalemated by inherent uncertainties of fluctuations 
called cosmic variance. Directional effects arc much more decisive, because 
they confront a symmetry, de Oliveira-Costa, ct al (2004) constructed an 
axial statistic for which they found modest statistical significance in multi- 
poles for I — 2,?> being rather well aligned. The fact that the dipole (/ = 1) 
also ahgns very closely with the multipoles 1 — 2,?, was later highlighted 
by Ralston and Jain (2004) and Schwarz et al (2004). When the dipole is 
interpreted to be due to our proper motion, it has sometimes been excluded 
as having "no cosmological significance." However there are many physi- 
cal mechanisms which can correlate CMB observations with galactic motion. 
There are good reasons not to exclude it, and the correlation of all three mul- 
tipoles significantly aligned with the constellation Virgo is quite inconsistent 
with chance. 

Subsequently there have been a large number of studies (Eriksen et al 
2004, Katz and Weeks 2004, Bielewicz et 0^2004, Bielewicz et aZ2005, Prunet 
et al 2005, Copi et al 2006, de Oliveira-Costa and Tegmark 2006, Wiaux et 
al 2006, Bcrnui et al 2006, Freeman et al 2006, Magueijo and Sorkin 2007, 
Bernui et al 2007, Copi et al 2007, Helling et al 2007, Land and Magueijo 
2007) which claim the CMB is not consistent with isotropy. Several physical 
explanations for the observed anisotropy have been put forward (Armendariz- 
Picon 2004, Moffat 2005, Gordon et al 2005, Vale 2005, Abramo et al 2006, 
Land and Magueijo 2006, Rakic et al 2006 Gumrukcuoglu et al 2006, Inoue 
and Silk 2006, Rodrigucs 2007, Naselsky et al 2007, CampaneUi et al 2007). 
It has also been suggested that the anisotropy may be due to foreground 
contamination (Slosar and Seljak 2004). Land and Magueijo (2005) find 
evidence that the detected anisotropy has positive mirror parity. Meanwhile 
some studies find no inconsistency ( Hajian et a/ 2004, Hajian and Souradccp 
2006, Donoghue and Donoghue 2005). There have also been several studies 
of the primordial perturbations (Koivisto and Mota 2006, Battye and Moss 
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2006, Armendariz-Picon 2006, Pereira et al 2007, Gumrukcuoglu et al 2007) 
and inflation (Hunt and Sarkar 2004, Buniy et al 2006) in an anisotropic 
universe. 

It is clearly necessary to explore new observables to determine whether 
anisotropy may be signs of physics beyond the standard paradigm. Here 
we develop new methods to test CMB data for isotropy. The new tests are 
possible because there exist many more invariant and vector- valued quantities 
than commonly examined. As conventional, the temperature distribution 
AT(n) is expanded in terms of the spherical harmonics, defining AT(n) = 
Yliim.'^irnYimin), whcrc n is a unit vector on the sky. The assumption of 
statistical isotropy is the statement that the ensemble average is given by 

< O'lmO'l'm' >= Ci5iii6mm'- (l) 

As a consequence all tensors formed from products of aim must be isotropic. 
We concentrate on the second rank tensors that can be formed, which include 
the isotropic power Ci 6ij and the power tensor Aij{l). They are defined by 

A{1) 6ij = 21 _^ I (hm^ij^mm'0'lm''i 
mm' 

^^ijO) = ^/^ _|_ ^\ '^Zm( JiJj )mm'0'lm'- (2) 

mm' 

Here Jj are the angular momentum operators in representation / for Cartesian 
indices i = 1...3. From Eq. [1] statistical isotropy predicts the ensemble 
averages 

< A>= Cu 
< >= CAj. (3) 

Giving attention to the power tensor represents a modest step towards 
examining the data more thoroughly. Consider the region of < / < 300, 
which includes the acclaimed peak in plots of /(/ + 1)C;. There are 299 data 
points from 2 < / < 300 but the information makes a smooth curve, one 
that can be fit with 3 or 4 parameters. The existence of an acoustic peak 
is useful but only the beginning of tests. Each aim has 2/ + 1 components, 
and the remaining 90,298 numbers, measured at great effort and considerable 
public expense, are never examined with the power spectrum. One reason 
comes from early tradition when data was poor. There also appears to be a 
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widespread but false belief that the power spectrum is the only rotationally 
invariant quantity. 

The approach of Copi et al (2007) consists of factoring the spherical har- 
monics expansion into products of vector-valued terms. All spherical har- 
monics can be developed from traceless symmetric polynomials of vector 
products. There are 21 + 1 real freedoms in multipoles of given I. Each mu- 
tipole can be factored into / unit vectors, each with 2 freedoms, making 21 
terms (Weeks 2004, Katz and Weeks 2004). Finally there is one overall power 
coefficient. The analysis of 50 vectors for / = 50 (say) is inherently compli- 
cated, and complicated further by discovery that Maxwell-multipole vectors 
of isotropic data are kinematically correlated to "repel" one another. Dennis 
(2005) discusses this and even solves the 2-point correlation analytically. 

Factoring into products of vectors is natural for tensors of rank 2. When 
we consider rank-2 tensors, we realize they have more than one invariant 
eigenvalue, so that the CMB power cannot even describe the quadrupole 
component. Tensors of rank higher than 2 have many types of invariants. 
One complication in finding them is the absence of a unique canonical form 
for high-rank tensors (da Silva and Lim 2006). Another technical obstacle 
is to obtain the invariants using quadratic functions of the data, which are 
statistically robust. 

The technical problems are solved here with a procedure exploiting an 
invariant canonical form where unique and preferred vectors are factored from 
each multipole. The fact this can be done is both obvious from symmetry and 
systematic via Clebsch-Gordan series. It is a particular linear transformation 
of the aim that leads to the power tensor Aij as quadratic form. When our 
procedure is applied to WMAP data there are numerous new invariants and 
several new signals of anisotropy. 

We had two distinct motivations in beginning our study: 

• Isotropy is a well-defined statistical data question that does not need 

any physical hypothesis to be motivated. Physical hypotheses are best 
characterized and tested through their symmetries. No free parameters 
are used in expressing isotropy, which can be confirmed or ruled out 
independent of the parameterization of model predictions. All of the 
new tests in this paper are also parameter-free. Isotropy of the CMB is 
also a separate issue from isotropy of other observables in the universe. 
This is because the CMB contains both information on formation in 
interaction with matter, and information on the propagation of light 
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since the era of decoupling. By separating tests of isotropy of the CMB 
we maintain an orderly process of testing larger issues in cosmology. 

• Anisotropy of the CMB may be a signal of physics beyond electro- 
dynamics and gravity. The interaction of light with an anisotropic 
axion-like condensate, or its conversion into weakly interacting parti- 
cles, may give anisotropic signatures without disturbing the large scale 
distribution of matter. 

Interestingly, there is a long history of "Virgo alignment" of electro- 
magnetic propagation effects which cannot be explained by any known 
media of conventional physics. Propagation of radiation in the sub- 
GHz region of radio telescopes, the 20-90 GHz region of CMB studies, 
and the optical region all reveal anisotropic effects. Many years ago 
Birch (1982) observed that the offsets of sub-GHz linear polarizations 
relative to radio galaxy axes were not distributed isotropically. The 
statistical significance of Birch's signal was confirmed by Kendall and 
Young (1984) and a statistic used by Bietenholz and Kronberg (1984), 
but generally rejected nonetheless. The redshift-dependent signal found 
by Nodland and Ralston (1997) was also generally dismissed. Rejection 
was based on not finding correlation in statistics believed to test for the 
same signals. Yet one must be sure the same tests were really made. 
Paradoxes were finally resolved by classification under parity. Statis- 
tics finding no signal test for even parity distribution with no sense 
of twist, while those statistics of odd parity consistently show signals 
of high statistical significance. Subsequent analysis (Jain and Ralston 
1999, Jain and Sarala, 2006) has found a robust axial relationship of 
odd parity character and aligned with Virgo. 

Until recently these explorations of anisotropy were conducted in the ra- 
dio regime, which might indicate a conventional explanation in plasma 
electrodynamics. The discovery of optical frequency polarization cor- 
relations (Hutsemekers (1998), Hutsemekers and Lamy (2001)) from 
cosmologically distant QSO's is extraordinary. The optical correlations 
again generate an axis aligned with Virgo (Ralston and Jain 2004). The 
coincidence of 5 closely aligned axes coming from 3 distinct frequency 
domains makes an exceptional mystery that merits further study. 

In this paper we make no attempt at physical explanation of anisotropy. 
Our new methods provide two independent new roads for data analysis. Or- 
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der by order for each / > 2 there exist three rotationally invariant eigen- 
values of Aij{l). The sum of three eigenvalues reproduce the usual power 
Ci. Meanwhile two independent combinations are new invariants. For the 
region of 2 < / < 50, say, we have 98 new invariants never examined be- 
fore. The isotropic CMB model predicts that all three eigenvalues should 
be degenerate and equal to C//3, up to fluctuations. We test this prediction 
by introducing the power entropy of Aij. Inspired by quantum statistical 
mechanics, the power entropy is a stable and bounded measure of the ran- 
domness of eigenvalues of matrix Aij. The entropy of the data turns out not 
to be consistent with the conventional CMB prediction. Next we examine the 
eigenvectors of Aij. The eigenvectors are independent covariant quantities. 
Under the isotropic prediction they should be randomly oriented. For the 
range of 2 < / < 50 there are 98 independent eigenvectors, and 49 "principal" 
eigenvectors associated with the largest eigenvalues. The largest eigenvalue 
makes the largest contribution to the total power and gives an objective rea- 
son to be singled out. Statistical measures of the ahgnment of eigenvectors 
do not happen to support the random isotropic expectations. Instead 6 axes 
obtained from the range 2 < / < 50 are found to be well aligned with the 
dipole or Virgo axis, each case having independent probability of less than 
5%. 

1.1 Outline of the Paper 

Transformation properties and signal processing of aim are encoded by linear 
operations. To take advantage we first introduce ip^, which is the unique 
linear transformation of aim that divides each Z-multiplet into direct products 
of irreducible representations with vectors. We relate ip^^^ to the power tensor, 
and study the distribution of invariants in the ILC map. 

Complementary to these studies are comparisons with "masked-sky" power 
spectra, evaluated with a window-function blocking out a particular angular 
region. Our linear transformations are made to mesh perfectly with the lin- 
ear operations of window functions and subtractions. In studies we report, 
we sometimes masked out the plane of the galaxy, and divided the remainder 
into two "top" and "bottom" regions. We study statistics of the regions as 
given, and also the statistics of full-sky maps they predict under standard 
procedures. There are several reasons for the masked sky studies. First, we 
can both study top/bottom asymmetries as a simple probe of isotropy in- 
cluding the galactic plane. We also study "cone-masks" which retain only a 
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region subtended by cones at high galactic latitudes, reducing the importance 
of galactic plane contamination. Finally wc make simulations of galactic con- 
tamination and show that our new methods readily detect contamination and 
find its angular location. The ILC data maps turn out to show significant 
anisotropy that cannot be attributed to galactic foregrounds. 

In the following section we describe the new methods. Applications to 
data follow. Not surprisingly, the power tensor method is generally more 
sensitive to anisotropy than power spectrum methods. Yet power spectrum 
methods find similar signals of anisotropy. This is significant because the two 
approaches are statistically independent - nothing about the power spectrum 
is used in the new approach. Almost all of our studies develop new informa- 
tion that by construction is independent of the usual CMB power 

The reliability of CMB maps is not under our control and certainly open 
to debate. The observers of the ILC team present maps claimed to have 
systematic errors equal to or smaller than cosmic variance. For this reason 
we don't discuss systematic errors. Our results are quantified using statistical 
errors based on Monte Carlo simulations with random realizations in place 
of the aim- 

2 Covariant Frames and Invariant Singular 
Values 

In this Section we use Dirac notation with 

aim ^ {l,m\ a{l)) . 

(There should be no confusion with ensemble averages.) The states \lm) are 
eigenvectors of angular momentum and with eigenvalues m, I {I + 1). 
The temperature measurements are real valued, which produces a constraint 
in the usual basis convention: 

aim = (-1)" a^_m 

Although the phase conventions standardize the aim they have no effect on 
anything we report. 

Besides simplicity, the power tensor is well motivated on rotational sym- 
metry grounds. Each multipole aim ~^ \a{l)) represents a geometrical object 
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visualized as a complicated beach-ball. "Holding by hand" the object for 
inspection, make a small rotation by angle 9, under which 

|a(/)) ^ |a(/))'= |a(0) + \Sa{l)) ; 

\6a{l)) = -iO- J\a{l)). 

Under what axes does the rotation make the most difference ? Compute the 
Hessian of the change, 

{5a{l) \ 5a{l)) = {a{l)\ f \a{l)) 



By Rayleigh-Ritz variation, the maximum rotation is developed along the 
eigenvectors of A^^ , which are natural principal axes of the aim ^ objects 
order by order in /. 

It is convenient to define a linear map or "wave function" ip^i^) (Ralston 
and Jain 2004) we call "vector factorization:" 

^'Jl) = m\j'\a{l)). (4) 

The purpose of this transformation is to extract (algebraically divide out) 

a vector (spin -1) quantity from each representation of spin—/, so each can 
be covariantly compared across different /. Inserting a complete set of states 
gives 

1 ' 

i^'rail) = -^=_ 5^ (/,m|J^|/,m')(/,m'|a(/)) 
I 

= ^mm'<^lm' (5) 

m'=—l 

where 

^im' - -^^=^^{1, ^\J% m') (6) 

Under rotations each index of ip^ rotates by its representation R{j), namely 

i^t^i^'^-R'^'mmmWl (7) 
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The transformation a ^ ip is invertible. The inverse relation can be written 

— \ ^ -pmm' i k 
'™ ~ k rm ) 

km' 

where 



-pk 

mm' ■ 

Upon rotating ip by the rule of Eq. [71 then 

a«m a'lm = T/T™ -R'''' (l)-Rmm'(0V'm5 

The standard power spectrum invariant is one of the invariants produced 
from ip^: it is 

mfc m 

Mathematically ip represents the multiplet a as a unique sum of outer 
products of vectors e" (spin-1 representation) with representations of 
spin—/. The terms are obtained by the singular value decomposition (SVD), 

3 

a=l 

Here and in subsequent equations we suppress label / when it is obvious. The 
Appendix discusses factorization in general terms of a Clebsh-Gordan series. 

The singular values A" are invariants under rotations on the indices k 
and m. (Indeed they are invariants under the even higher symmetry of 
independent 5*0(3) x 5*0(3) rotations.) The various factors are constructed 
by diagonalizing the 3x3 Hermitian matrices 

m 

= E er(A")'er; (9) 

a 

(V^V)mm'(0=E^-(^)^-'(0, 

k 



Since they are the eigenvectors of Hermitian matrices the and are 
generally orthogonal, except for exceptional degeneracies, and normalized to 

k 
m 

Since (A")^ are the eigenvalues of Hermitian matrices they are real and pos- 
itive, with A" > defining the sign convention for e". Since aim are real, we 
also have e" real valued. 

The set of three orthogonal e" defines the "frame vectors", which make 
a preferred frame for the vector components of ip. In that preferred frame 
and the frame of three m", the matrix ip^ is diagonal with three invariant 
singular values A°. We will call Aa6ai3 the ^^SV matrix", SV standing for 
singular values. 

2.1 Isotropy 

The relation of to the power tensor Aij is simple algebra. We have 
A^^{l) = j^^tr{.r\a{l)){a{l)\jn, 

where tr is the trace. 
Then 

m 

a 

which is just Eq. [9l Isotropy holds that the eigenvectors e° must be dis- 
tributed isotropically on the sky and that the eigenvalues A"" are random 
variables concentrated at a/G/S. In terms of the V'm factors, Eq. [T] predicts 

< >= Y^''5mm'6'''. (10) 

Eq. [To] can be contrasted with the correlations of Maxwell-multipole 
vectors developed by dividing multipoles into products of vectors. Those 
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vectors have not been organized into irreducible representations nor classi- 
fied in invariant terms of importance, and they do not have an uncorrelated 
distribution (Dennis 2005). An advantage of using our ip^ representations is 
that isotropy transforms to isotropy without annoying Jacobian factors. 

Here is how to interpret the factors. The highest possible degree of 
anisotropy produces one singular value equal to the total power, and two 
others that vanish. The corresponding aim components can be written as the 
product of one vector e'-^^ and one multiplet ■u'^^\ We will call this special sit- 
uation a "pure state." This has a very simple realization for the quadrupole 
case 1 — 2. Any quadrupole is equivalent to a symmetric 3x3 matrix made 
from sums of products of 3-vectors. A generic symmetric matrix has 3 real 
and unequal eigenvalues. If two eigenvalues vanish then the matrix is the 
outer product of a unique vector, a pure state. Conversely, if all eigenval- 
ues are degenerate, the matrix is a (trace subtracted) multiple of the unit 
matrix, equivalent to the sum of outer products of three frame vectors by 
completeness, 1 = |e") (e"^| , which is the isotropic prediction. 

2.2 Tensor Power Entropy 

The isotropy hypothesis that the SV matrix should be \fC^ times the unit 
matrix can be tested in several different invariant ways. 

The information entropy or simply power entropy of the power tensor 
comes by recognizing that p''^ = {jpilj'^Y^ is proportional to a density matrix 
on 3-space. The proportionality constant is the overall power. Removing it 
produces a normalized form 



Von Neumann (1932) found the entropy S of normalized Hermitian matrix 
p to be 





(11) 



s 




(12) 



^(A")2%((A«)2), 



a 







< 



S < log{3). 



Here (A")^ are normalized to sum to one, again removing the overall power 
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from discussion: 

This normalizes tr{p) = 1. 

The entropy is unique in being invariant, extensive, positive, additive for 
independent subsystems, and zero for pure states. A density matrix with 
no information is a multiple of the unit matrix, and has entropy equalling 
the log of its dimension. Very simply Siso log{3) is the isotropic CMB 
prediction. Low entropy compared to isotropy is a measure of concentration 
of power along one or another eigenvector of p. In the present context pure 
states have all the power in one singular value, define one single directional 
eigenvector, and have power entropy Spure —* 0. 



2.3 Power Alignments Across /-Classes 

Our methods allow us to explore the alignment of power tensors between 
different /-classes. 

One of the most interesting tests concerns the "principal axis e(/)", which 
means the eigenvector with the largest eigenvalue for each /. From Eq. [TO] the 
set of all principal frame vectors should be a symmetrical ball in the isotropic 
prediction. If the data is anisotropic a bundle of vectors may lie along an 
axis or in a preferred plane. It is important to remember that the sign of 
eigenvectors is meaningless, and determined by algorithms assigning signs 
in computer code. Thus the "average" eigenvector is not a good statistic. 
Tensor products are the natural probe of a collection. For statistical studies 
we construct a matrix X defined by 

1=2 

The eigenvalues of X are a probe of the shape of the bundle collected from 
2 < / < ^max- We probe the isotropy of the collection with the directional 
entropy Sx computed using Eq. [12] and X px- The directional entropy 
does not use the singular values of the CMB data, and is independent of the 
power entropy. Confirmation of isotropy comes if Sx ~ log{?)) up to random 
fluctuations. A signal of anisotropy would be an unusually low value of Sx 
compared to log (3). 
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2.3.1 Traceless Power Tensor 



We also compare alignments using a statistic that includes the weighting by 
singular values. For this purpose we define the traceless power tensor i?*-', 

B'^{1) =A'^{1) - -tr{A{l))6'^. 
3 

The eigenvectors of B are the same as A, but the value of the total power 
has been removed. To compare two angular momenta classes we examine the 
correlation 

^ tr(B(mn) ,15, 

y'tr{B(l)B{l))y'tr(B{l')B(l')) 
In the isotropic uncorrelated model Y{1, I') should be proportional to 6ii/. 



3 Full Sky Studies 

The WMAP-ILC team developed a foreground cleaned temperature anisotropy 
map using a method described as Internal Linear Combination (ILC). The 
process combines data from 5 bands of frequencies 23, 33, 41, 61, and 91 
GHz. While there also exist several other cleaning procedures that are inter- 
esting to compare (Tegmark et a/ 2003, Saha et a/ 2006, Eriksen et a/ 2007), 
our main focus lies on the ILC map. 

According to Hinshaw et al (2007) the ILC map is known to become 
dominated by statistical instrumental noise for / > 400. Errors on the aim 
come from several sources. Systematic errors may occur from inappropri- 
ate foreground subtractions. Ambiguities in combining frequency bands also 
contribute. The range 2 < / < 50 is considered to be reliable under state- 
ments that statistical and systematic errors lie within ranges typical of cosmic 
variance. We restrict our studies to this range. 

Fig. d] shows the normalized eigenvalues A", with a dashed line for the 
isotropic prediction Aq, = Cil\f?>. The normalization is (A")^ = 1. The 
power entropy Sfua-sky{l) mode by mode in / for the ILC map is shown in 
Fig. [21 By Monte Carlo simulations we find that the points 1= 6, 16, 17, 30, 
34, 40 appear to be statistically unlikely. 

We generated random aim and derived the power entropy using 10,000 
realizations of CMB maps. In our procedure real aio values were drawn from 
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Figure 1: Normalized singular values versus / for the ILC full sky data set. For 
each / there are 3 values A"^ normalized so the sum of their squares is unity. The 
dashed line shows the isotropic prediction A" = l/\/3- 
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Figure 2: Power entropy S{1) of the ILC full sky singular values versus mode 
number I. Low power entropy indicates orderly behavior not typical of isotropic 
data. The dashed line corresponds to the maximum values of the entropy S = 
log{2,). 
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Gaussian distributions with zero mean and unit norm. Complex for 
m = 1, were created by adding real and z-times real numbers from the 
same distribution and dividing by \pl. As a consistency check, P-values were 
computed independently by different members of the group, both including 
the Ci values, and omitting them, to verify the method and that the usual 
power statistic indeed drops out. Finally we computed P-values representing 
the frequency that power entropy in a random sample is less than that seen 
in the map. In Fig. [3] we show P-values of power entropy realized in the ILC 
maps versus mode number /. There are 6 unusually low power entropies of 
5% or smaller found at / =6, 16, 17, 30, 34 and 40, with P-values of 0.040, 
0.032, 0.041, 0.018, 0.045, 0.024 respectively. 
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Figure 3: Zo5fio(-P)-values of the full-sky ILC map power entropy shown in Fig. O 
P-values are estimated from 10,000 realizations of random isotropic CMB maps. 
The dashed horizontal line shows P = 5%. 

It is interesting that the power entropy does not single out / = 2 or 3 as 
particularly significant. The power entropy itself does not have any informa- 
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tion about the directional alignment of eigenvectors, which is independent 
and will be studied separately. 

Note there are many independent instances of P-values less than the ex- 
pected average of 1/2. What is the probability to see so many small P- 
values? Since the P-values already take into account the distributions, this 
basic question has an analytic answer. 

There are 6 power entropies with P values of p = 0.045 or less. The 
distribution to get k independent P-values, found less than or equal to p, 
among n independent trials is /( A; | n, p) = p'^(l —p)"'~^n\/{n — k)\k\ This is 
the binomial distribution, which comes from counting k instances of "passes" 
when considering "pass-fail" questions, while accounting for n\/{{n — k)\k\) 
independent ways to get the outcome "pass." The total number of indepen- 
dent trials with /^ax = 50 are n = 48. The total probability of the data given 
random chance i^j 

/( 6 events | n = 48, p = 0.045) = 0.015 (16) 

This clearly shows that the signal of anisotropy is present over a much larger 
/ range in the CMB data in comparison to what is commonly believed in the 
literature ( Katz and Weeks 2004, Bielewicz et al 2004, Bielewicz et al 2005, 
Prunet et al 2005, Copi et al 2006, de Oliveira-Costa and Tegmark 2006, 
Bernui et al 2006, Freeman et al 2006, Magueijo and Sorkin 2007, Bernui et 
al 2007, Copi et al 2007, HeUing et al 2007, Land and Magueijo 2007). 

It is interesting to throw out unlikely events one by one. The binomial 
probability to find 6, 5, 4... small P-values sorted from 0.045 to smaller levels 
is 1.5 X 10"^, 3.1 X 10"^, 7.5 x 10~^,... respectively. We also conducted a 
Monte Carlo simulation to verify the distribution is as stated. 

A common criterion to test hypotheses poses confidence level break-points 
of 5%. The distribution of the invariant power entropy is not consistent with 
the isotropic model. 

3.1 Alignments 

The alignment of eigenvectors in the ILC Map at low / may be a signal 
of a fundamental anisotropy or could also be caused by galactic foreground 

^In discrete data analysis each particular possibility is enumerated, as in throwing 
dice, so we put no emphasis on the cumulative distribution. The cumulative binomial 
probability to see 6 or more cases is slightly larger, /(fc > 6 events \ n — 48, p — 0.045) ~ 
0.020. 
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Figure 4: /o5fio(P) -values of largest power eigenvectors to align with the 
quadrupole versus /. The dashed line shows the 5% level of significance. 
There are 6 points at the level of P < 5% or less. 



contamination. The history of alignment with Virgo suggests we compare 
the principal axes of the power tensors /-by-/ with the quadrupole axis. 

Fig. m shows /o(?io(-P)-values in an isotropic distribution for power prin- 
cipal eigenvectors to align with the quadrupole. There are 6 axes aligned 
with a significance close to 5% level or less; they are multipoles with / = 
3, 9, 16, 21, 40, 43. The binomial probability of seeing 6 such events is 
2.3%. Hence the alignment with the quadrupole is seen over a relatively 
large range of / values. If coincidences are thrown out one by one, the es- 
timated probability to find (6, 5, 4...) small P-values seen in the data is 
2.3 X 10~^, 1.7 X 10~^, 1.4 X 10"^.... The distribution of the principal axes 
for many different / are not at consistent with the isotropic proposal. The 
Cartesian components of the well-aligned principal axes are listed in Table 

m 

It may appear arbitrary to compare the axes to the / = 2 case. It is 
motivated by the literature and our previous history of work, discussed in 
the Introduction. For completeness P — values among all the independent 
pairs are listed in Table [2l It shows that choosing the / = 2 case causes no 
special bias, because once the axes are well-aligned any one of them could 
be used for comparison. Correlation with the dipole, which is not strictly 
part of the ILC map, is also included for reference in Tabled! The statistical 
significance against isotropy is high whether or not one includes the dipole. 
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/ 


axis h { X, y, z ) 


2 


(-0.209, -0.302,0.930) 


3 


(-0.251, -0.385, 0.888) 


9 


(-0.224, -0.0974, 0.970) 


16 


(-0.157, -0.175, 0.972) 


21 


(-0.139, -0.536, 0.832) 


40 


(-0.0928,0.211,0.973) 


43 


(-0.217, -0.132, 0.967) 


1 


(-0.0616, -0.664, 0.745) 


nx 


(-0.435, 0.134, 0.890) 



Table 1: Galactic Cartesian coordinates of principle axes ("large eigen- 
vectors") n that are extraordinarily aligned. Criteria for selection are 
P — values < 5% relative to the / = 2 principal axis. The / = 1 or dipole 
axis is included for completeness. The axis fix is the principal eigenvector of 
the directional entropy of the set 2 < / < 50. 





/ = 2 


3 


9 


16 


21 


40 


43 


r = 2 
















3 


0.005 














9 


0.022 


0.045 












16 


0.010 


0.030 


0.005 










21 


0.035 


0.019 


0.109 


0.075 








40 


0.051 


0.078 


0.057 


0.032 


0.090 






43 


0.015 


0.036 


0.0006 


0.003 


0.094 


0.051 




1 


0.094 


0.067 


0.199 


0.150 


0.015 


0.141 


0.178 



Table 2: P— values of coincidence between independent pairs of principal axes 
labeled by I, I' shown in Table [TJ Correlation with the I' = 1 axis, which is not 
strictly part of the ILC map, is included for completeness. 
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3.1.1 Alignment Entropy 

In performing tliis researcli we first tested for clustering using the directional 
entropy Sx of matrix X. Comparison with a Monte Carlo simulation using 
the same range 2 < Z < 50 does not show anything extraordinary, yielding 
a P-value of 69.0 %. The directional entropy Sx did not signal significant 
clustering. Yet the directional entropy is rather insensitive statistic which 
misses the pattern in Table [2l 

Curiously, however, the principal axis of X is 

fix = (-0.435, 0.134, 0.890). 

The angular galactic coordinates (/ = 162.88, b = 62.91) are clearly not in 
the galactic plane, but well aligned with Virgo. Recall that the principal vec- 
tors of the quadrupole and octupole in the ILC Map have galactic Cartesian 
components 

n(/ = 2) = (-0.209, -0.302, 0.930) 

n{l = 3) = ( -0.251, -0.386, 0.888 ) (17) 

respectively. These two axes are very closely aligned with one another, as so 
much noted. It is remarkable that they are aligned with the principal axis of 
X coming from the whole ensemble. 

Our first studies also used a Mollweide projection to examine visually for 
clustering of axes (Fig. [6]). Since we are plotting an axis, and not a vector, 
we plot points only on one-half of the sphere. The ILC map does not visu- 
ally show any striking clustering among different eigenvectors in Mollweide 
projection. 

However the Mollweide projection seriously distorts the region near the 
galactic poles. Fig. [7| shows a stereographic projection, a method that is 
free of distortion near the poles. Here the clustering of axes is clearly visible, 
explaining that the Mollweide projection has sent all the interesting points 
out of view. 

3.1.2 The Region 1 < / < 11 

The method of Copi et al (2007) finds highly significant signals for the region 
1 < / < 11. It is interesting to evaluate our results for this range. 

There are /our well-aligned principal vectors at Z = 1, 2, 3, 9. Removing 
one as trivial, the probability to find 3 p- values of less than 5% probability in 
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Figure 5: Directional entropy histograms of Sx for the multipole range 2 < I < 
50. The value of directional entropy of the ILC map are also shown. Histograms 
were generated using 10,000 randomly generated CMB data sets. 




Figure 6: Mollweide projection plot in galactic coordinates of principal axes for 
the ILC map over the range 2 < / < 50. No clustering is visible because the cluster 
of axes near the galactic pole is distorted. 
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Figure 7: Clustering of principal axes visible in stereographic projection from 
the North galactic pole. Highly correlated axes (Table [Hare circled, with the 
quadrupole reference point (large circle). Also shown are the dipole 

(diamond), the radio polarization offsets alignment axis (star), and the QSO 
optical polarization alignment axis (triangle), as discussed in the text. The 
inset box shows an enlargement and 5° separation scale in the cluster region. 
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a random sample in 11 trials is about 1%. One unusually low power entropy 
value exists in the range at / = 6. 

What about the directional entropy ? When Sx is evaluated for the range 
2 < / < 11, random samples show Sx equal or lower than the ILC Map only 
440 times out of 10,000, a P-value of AA%. This is another indication of 
anisotropy. The P-values for different upper limits /upper are shown in Fig. [HI 
We see that most of the P-values in this range are smaller than 5 %. Finally 
the principal axis of X is 

nx(2 < / < 11) = (-0.531, 0.117, 0.839). 

It is well-aligned with Virgo. The signals from 1 < / < 11 evidently come 
from the same source as the signals from the entire range 2 < / < 50. 

3.1.3 The Y Matrix 

The y-matrix gives another method to examine visually the angular corre- 
lations. Fig. [HI shows contour maps (color online) of the y(/, /') correlation. 
The correlation in ideal uncorrelated data should show a diagonal line up 
to fluctuations. Features of Fig. [9] confirm the other studies of unusual cor- 
relations in nominally uncorrelated random data. This study complements 
the study with matrix X, which uses principal axes and does not use the 
information contained in the singular values of ip. 

3.2 Foreground Contaminated Maps 

Here we study effects of simulated foreground contamination on the power 
tensor and entropy statistics. 

A contour plot shows the input of the simulation in FigJTUl top panel. The 
map is made by adding a random CMB map with a simulated synchrotron 
foreground map (Giardino et al 2002) generated at 23 GHz frequency and 
normalized to 2% of the actual foreground strength. Contamination is dom- 
inantly in the galactic plane using the procedure of . 

Our procedure readily senses the contamination. FigUHl bottom panel, 
shows the spatial distribution of largest eigenvectors e/ for the range of mul- 
tipole moments 2 < I < 50. The vectors lie dominantly in the galactic plane, 
providing a clear signal of alignment caused by galactic contamination. 

Calculation of the entropy of the largest eigenvectors Sxilmax = 50) = 
0.995 also shows significant anisotropy. The P value for this to occur in a 



23 



10 20 30 40 50 



10 20 30 40 50 



-1 



Figure 9: Contour maps of the correlation Y{1, I') of traceless tensor power 
between different angular momentum indices /. Left: a simple density plot. 
Right: Contours interpolated with trends more readily visible. A scale (color 
online) shows the color code of correlation values. Features confirm the other 
studies of unusual correlations in supposedly uncorrelated random CMB 
data. 



random isotropic sky is P = 0.01 %. The principal axis, or largest eigenvector 
of X is ri = (—0.233, —0.963, —0.132), which nicely corresponds to (Z = 
76.4, b = 7.56), lying close to the galactic plane. 

These studies show that our procedure detects galactic contamination in 
a simple and statistically reliable way. They also indicate that the Virgo 
alignment observed in the preferred directions of the ILC data cannot be 
attributed to in-plane galactic contamination. 

4 Masked Sky Studies 

As mentioned in the introduction, our use of linear transformations allows 
the effects of window functions of masks to be handled by conventional and 
transparent means. In this Section we explore the effects of masks applied to 
different hemispheres of the galaxy. Although masks distort power spectra, 
using symmetrical masking creates no biases. 

We use four different masks. The Northern and Southern hemispheres in 
galactic coordinates are retained in Mask 1 (Ml(ri)) and Mask 2 (M2(n)), 
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Figure 10: Top panel: Simulated CMB + (2%) synchrotron contaminated fore- 
ground map described in the text. The simulated foreground contamination is 2% 
of the actual foreground strength. Bottom panel: Principal axes for 2 < / < 50 
made from the map of the top panel show alignment in the plane of the galaxy. 
Mollweide projection in galactic coordinates. 
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respectively. Cone Mask 1 (CMl(n)) and Cone Mask 2 (CM2(n)) respec- 
tively retain the polar angle regions < ^ < 30° and 150° < 9 < 180°. These 
are typical values chosen to minimize galactic contamination. We show the 
ILC Map cut by Ml and CMl in Fig. [H 

4.1 Isotropic Maps Estimated from Masked Maps 

We explore the power estimated for a full-sky map, as commonly developed 
from the power spectrum obtained from unmasked regions. A window func- 
tion W{h) is defined so that 

AT""^^'(^) = W{n)AT{n). 

This relation is linear and invertible, and already in its diagonal form: 

Ar(fi) = ^T'^^'^{n)/W{h) (18) 

Inversion is unstable ("ill-conditioned") for small W{n), the region where 
information does not exist, a fact that naturally limits to information that 
can reliably be obtained. 

We used the method of Hivon et al, as we briefly review. Let aim{mask) 
be the amplitude obtained from AT"^""'^^ (n) . The transformation of Eq. [18] 
can be carried out in angular momentum coordinates, with the general form 

aim{mask) = yv'^^,ai'm', 
|a(/)(masA;)) = W"' |a(/')), (19) 

where W is the transformed kernel. In general all components of aimifnask), 
and not just the power spectrum are needed. The power spectrum is related 
by 

{a{l){mask) \ a{l){mask)) = {a{l")\ W"'*W'' W)) . (20) 

I 

An estimate of the power spectrum from inverting Eq. [20] will be linear but 
generally couple all multipoles in complicated ways. 

To construct an inversion Hivon et al make the assumption of isotropy. 
In Eq. [20] replace 

(m I a(/)) (a(/") | m") ^ Q5""<5w. (21) 
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Figure 11: The WMAP-ILC map after applying masks Afl(n) and CMl(n). 
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Assuming this applies to ensemble averages, the linear relation between power 
measures is: 



< Ci{mask) >= J] Mu, < Cv >, (22) 



where M/^/j is computed with Clebsch (3j) coefficients 

M.. = ^T'p'3+m('- '^ X (23) 

Here Wi is the power spectrum of the mask function. Beam and pixel factors 
are readily incorporated in Eq. [221 while for our purposes they can be ab- 
sorbed into the symbols. Since Eq. [22] is linear, the kernel M can be inverted 
to solve for estimated C; of an ideal full sky in terms of the observable power 
C/(masA;). 

In general inversion remains ill-conditioned, and care must be used in 
interpretation. Hivon et al choose to work with binned power spectra defined 
by a running average: 

Cb = Yl PbiCiimask). (24) 



The "binning operator" Phi we use averages 3 multipoles, with 

I'' 

2tt 



'■max 7/7 _|_ 1 N 



. Binning may ameliorate correlations caused by windowing, while generally 
creating its own correlations between different Ci,. The estimated full sky 
power spectrum Cb' is 

Cb' = E KylCb] 
b' 

Kbb' = (P • M • Q)bb' 

The "Inverse binning operator" Qib is given by Qib = j^f^ and is nonzero 
only if 2 < l^^^ < I < I'^ax- This summarizes sophisticated data processing 
procedures applied to the BOOMERANG experiment and many others. 
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We tested our code for full sky power spectrum estimation using Monte- 
Carlo simulation. We generated 100 random CMB maps, masked them by 
each of the four masks, and obtained partial sky power spectra. We binned 
the average masked sky power coefficients to decrease dimensions from an 
original size of 1024 x 1024 to a reduced size of 341 x 341. Finally we produced 
an estimated full sky power. Figure [12] shows the results of the simulations. 
The estimated full sky power spectrum matches quite well with the average 
of the input full sky binned power spectra for each of the four masks. 

We note that the binning and isotropic assumptions of the full sky es- 
timation method naturally dilute statistical measures. We did not pursued 
the systematic errors. We simply applied our procedures and studied the 
outcomes. 

4.2 Features of Estimated Full Sky Spectra 

In Figure [TSl we show the estimated full sky power spectra obtained from ILC 
using the four masks Ml(n), M2(n) and CMl(n), CM2{n). We observe a 
definite asymmetry in the estimated full-sky power spectrum between the 
Northern and Southern hemisphere in several multipole ranges. 

It is interesting that Ml(n), M2(n) and CMl(n), CM2{h) find similar 
excess in the same multipole ranges. This shows that the power anisotropy 
observed cannot be explained away by contamination of the galactic plane, 
but also exists at high galactic latitude. To quantify the power differential we 
sum the difference in power over all bins. The probability of randomly getting 
an equal or larger difference is less than 1 % for the range 2 < / < 11 and 
about 9 % for the range 2 < / < 50. These figures are given retrospectively 
after scanning the figures, and one may feel they represent a search. Given 
the unknown systematic errors we will not bother estimating significance. It 
is nevertheless interesting that a strictly power-based search, which includes 
averaging features that tend to wash out signals, can still potentially detect 
anisotropy. 

Just for interest we also include the estimated full sky ILC power from 
an unhinned inversion (Fig. [T^ . There are substantial differences from the 
binned evaluation. For example, negative power occurs in some cases. Since 
the unbinned procedure retains more information about the inversion than 
binning, it is a symptom of the instability of inversion, and not a mathemat- 
ical error. Use of binning helps to cure the problem for an ensemble of many 
maps, but it weakens the relation of any particular map to any particular 
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Figure 12: Comparison of estimated full sky power spectrum with input power 
spectrum simulation using only masked regions. Blue (red) online denote the 
Northern and Southern Hemisphere cases respectively. Top panel: Mask 1 and 
Mask 2. Bottom panel: Cone Mask 1 and Cone Mask 2. The spectra match 
the average of the spectra generating them (solid black line, "Theory") nearly to 
the upper limits of power developed in the maps, / ~ 300. The inset shows the 
difference Diff. 
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Figure 13: Estimated full sky binned power spectra, Northern and Southern 
galactic hemisphere, using Ml{h) and M2(n) (upper panel) and CMl{h) and 
CM2{h) (lower panel). 
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MultipoleJ 



Figure 14: The estimated full sky power map constructed from Mask 1 with- 
out going through the running average ( "binning" ) procedure that combines 
different Q. Note that negative power is sometimes predicted. 



data. Exploring such features gives an idea of the reliability of trying to 
improve instabilities of inversion by binning. 



5 Summary and Conclusions 

We have developed several new symmetry-based methods to test CMB data 
for isotropy. 

The wave function ipmi^) is a linear transformation of the amplitudes 
aim which expresses them as products of vectors and representations /. The 
power tensor Aij is quadratic in the wave functions and aim- Invariants of 
the power tensor are its eigenvalues, the singular values of ip. The statistical 
distribution of these invariants is a test of randomness and isotropy of the 
the power tensor / by I. The power entropy serves as a robust statistic. 
The power entropy is independent of the usual power, and also independent 
of the eigenvalues of Aij, which contain information on the orientation of 
multipoles. To assess multipole orientations, we collected the principal axis 
of Aij, which is the eigenvector with largest eigenvalue. 

Both the power entropy and the axial vector distributions show features 
that are not consistent with isotropy. Fig. [TBI summarizes some results of the 
ILC-full sky study. The figure repeats Fig. [T] highlighting the exceptional 
cases. Cases of very low power entropy are joined by a short thing line. 
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Figure 15: Normalized singular values versus I with exceptional cases highlighted. 
Thin short lines denote P-values of entropy Spower ^ 5%. Thick gray lines denote 
principal axes aligned with the quadrupole (Virgo) axis with P < 5%. There are 
12 independent cases of 5% or lower significance, which is unlikely in isotropic 
data at significance exceeding 99.9%. Significance is increased when the dipole is 
included. 



34 



Cases of very high ahgnment of the principal axis with the quadrupole axis 
are shown by a longer thick gray line. There are 12 independent cases hav- 
ing 5% probability or lower in an isotropic distribution, with 2 coincidences 
showing both exceptional power entropy and alignment. Ignoring the double 
coincidences, there are 98 independent trials in examining 2 < Z < 49 twice. 
The total probability to see the data's values in uncorrelated isotropic CMB 
maps is of order 10~^. 

Our study substantially extends previous studies that find low / multi- 
poles of current CMB maps are highly unlikely on statistical grounds. Unlike 
methods that produce a large number of vector factors for each I, our wave 
function method identifies a single invariant axis for each I, and 3 indepen- 
dent invariants for each /. The axial quantities are remarkably aligned with 
the direction of Virgo in rather systematic fashion for many /. CMB data 
show 7 total axes afigned over the range 1 < / < 50, with each case hav- 
ing independent probability at order 5% or smaller. The values of I— 6, 16, 
17, 30, 34, 40 have anomalously non-random power eigenvalue distributions 
giving exceptionally lower power entropy. Multipoles with / = 3, 9, 16, 21, 
40, 43 have anomalous alignment with the quadrupole. It is no longer possi- 
ble to restrict anomalous CMB statistics to the relative orientation or power 
seen in the dipole, quadrupole and octupole cases. Our study shows, for 
the first time, that the alignment of the multipoles with an axis pointing in 
the direction of Virgo is much more pervasive and not just confined to low I 
values. 

To explore the origin of these features, including possible foreground ef- 
fects, we repeated many calculations using sky-masked data sets. The data 
shows evidence for systematic differences between the Northern and South- 
ern regions of the sky. Masking out the galactic plane does not eliminate 
signals of anisotropy seen in the full sky studies. As consistency checks, 
the anisotropics of CMB plus simulated synchrotron emission contamination 
from the plane of the galaxy are detected by our methods in just the spatial 
regions where they are simulated. The observed anisotropies cannot readily 
be explained away by appealing to galactic foreground contamination. 
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6 Appendix 



Here we relate the factorization wave function ijj^ to covariant canonical 
forms of the multipoles aim- 

In the theory of angular momentum one conventionally decomposes prod- 
ucts of representations of dimensions ji, j2 into sums of states \JM), where 
\ji — j2\ < J < ji + J2- It is also possible to divide the large representation 
\JM) into products of smaller representations. Let = (JM | be the 
matrix elements of a state |^) transforming on representation J. We write 

where the Clebsh-Gordan coefficients are 

r^i-u.-2* ^ ^j^rn,j,m2 \ JM) . 

The Clebsch's dividing angular momentum J in products of ji = 1 and 
j'2 = J are simply the angular momentum generators J cast into the angular 
momentum basis. The relation of to ipmimi ^^^^ 1^) ~^ k(0)' ^^^^ — 
mi, converts one Cartesian index k into angular momentum m2, uses ji — 1 
and j2 — J — i- 

There exist a large number of decompositions for different choices of fac- 
torization. Each particular choice has a privileged status when viewed as the 
unique way to factor the multiplet into products of smaller spaces. There is 
no principle upon which to choose the smaller spaces, explaining the remark 
that no preferred canonical form exists for high-rank tensors. The relation 
of ip^ to the canonical form of Copi et al can be viewed as coming from two 
steps. First divide \a) (ji = 1) (ji = 1) ..(j2 = !)••• ® -{ji = 1), to 
develop the Maxwell multipoles. Next find a unique preferred vector which 
multiplies a tensor of rank I made from traceless symmetric products of the 
other vectors. That vector corresponds to ■0m- Alternatively can sim- 
ply be written down (Ralston and Jain 2004) as the vector available from 
symmetry. 
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